PHYSICAL REWIEV E Sumbitted 
On Conditional Statistics in Scalar Turbulence: Theory vs. Experiment 
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We consider turbulent advection of a scalar field T(r), passive or active, and focus on the statistics 
of gradient fields conditioned on scalar differences AT(R) across a scale R. In particular we focus 
on two conditional averages (V 2 T\AT(R)) and (|VT| 2 1 AT(R)). We find exact relations between 
these averages, and with the help of the fusion rules we propose a general representation for these 
objects in terms of the probability density function P(AT, R) of AT(R). These results offer a new 
way to analyze experimental data that is presented in this paper. The main question that we ask 
is whether the conditional average (V 2 r|Ar(_R)) is linear in AT. We show that there exists a 
dimensionless parameter which governs the deviation from linearity. The data analysis indicates 
that this parameter is very small for passive scalar advection, and is generally a decreasing function 
of the Rayleigh number for the convection data. 



I. INTRODUCTION 

The equations of motion in fluid mechanics, be them 
for the velocity field u(r,t) or for a scalar field like the 
temperature T(r, t), contain interaction terms like it- Vtt 
or u ■ VT and dissipative terms like vV 2 u or kV 2 T, with 
v and k being the kinematic viscosity and the scalar dif- 
fusivity respectively. Accordingly, when one attempts 
to derive a theory of correlations functions of the field 
or of field differences across a length scale R one runs 
into mixed correlation functions of the Laplacian of the 
field with the field itself. At this point there are two 
fundamentally different approaches that we want to ex- 
pose first. We exemplify these approaches in the context 
of turbulent advection, but similar considerations apply 
also to Navier- Stokes turbulence. 

One fundamental strategy, which is the more usual one, 
is to consider the structure functions of the field differ- 
ences. Denoting AT(r, r') = T(r') — T(r), the structure 
functions S n (R) are defined by 



S n (R) = ([AT(r,r')} n ) , R=\R\ = \r'-r\ 



(1) 



where homogeneity and isotropy of the ensemble were 
assumed. Using the advection equation 



dT 



u-VT = kV 2 T 



one can derive the equation of motion for S n (R): 
dS n (R) 



Of 



+ D n (R) = J n {R) , 



D n (R) = 2n([u(r) ■ V AT(r)][AT(r, r')] n_1 > 
J n (R) =-2n K (V 2 T(r)[AT(r,r')]™- 1 ) . 



(2) 



(3) 

(4) 
(5) 



In the stationary state we have a "balance equation" 
D n (R) = J n (R). 

The obvious advantage of this approach is that it in- 
volves objects that depend on one coordinate R only. On 



the other hand the analysis of the balance equation re- 
quires a theory of the viscous term J n which involves a 
correlation between the Laplacian of the field and field 
differences. Even when R is within the inertial range, 
one cannot get rid of J n (R). The limit k — > does not 
help; J n has a finite value in this limit. 

A second fundamental approach which avoids this dif- 
ficulty 0] employs multi-point correlation function of 
field differences. Starting form the same field differences 
AT(ri,r[) one defines the correlation function 

^ n (n,ri, . . . r n ,r' n ) = (AT(n,ri) . . . AT(r n ,r' n )) , 

(6) 

which depends on n fields and 2n coordinates. The equa- 
tion of motion for T n looks superficially similar to (til): 



dt 



+ T> n (r 1 ,r' 1 , . ,.r n yj = J n (ri,r[, . ■■r n ,r' n ) , 



(7) 



where 



X>„(rx,ri; ...r n ,r' n ) = £<AT(n,ri) . . . [u{r' 3 ) ■ V^Tft) 

3=1 

-u(r j )-V j T(r j )]...AT(r n ,r' n )). (8) 

n 

Jnin,^; . ..r n ,r' n ) = kJ^ (V 2 + V 2 ,) (AT(n, r^) . . . 

3 = 1 

...AT(r j ,r' j )...AT(r n y n )) . (9) 

In the stationary state we again face a "generalized bal- 
ance equation" T> n — J n . In fact, there is a fundamental 
difference between the two approaches. In the present 
case one can analyze the generalized balance equation 
keeping all the separations within the inertial interval of 
scales. Then we can take the limit k — » and show jl] 
that the dissipative term J n vanishes in the limit. We re- 
main in this limit with a homogeneous equation T> n = 



1 



without any Laplacian terms. The advantage is that we 
have in principle a complete theory without the need of 
additional input. The obvious disadvantage of this ap- 
proach is that we have functions of many variables. Nev- 
ertheless, this approach turned out to be very useful in 
the context of Kraichnan's model for passive scalar advec- 
tion H , where the homogeneous equation can be turned 
into a linear partial differential equation for the correla- 
tion functions T n . But even in this simplest possible case 
the difficulty incurred by having functions of many vari- 
ables led to contradicting arguments about the relevant 
physical solution. 

It is thus obviously useful to find ways to analyze the 
simpler version in which we have one variable only, but a 
mixture of inertial and dissipative scales contributing to 
the correlation function J n . To proceed we need however 
additional information. One possible way of inputting 
this information is in the language of conditional aver- 
ages. To see this consider the mean value of the kV 2 T 
condition on AT: 

H{AT,R) ee n(v 2 T(r) AT(r,r')) . (10) 
Using this definition we rewrite J n {R) as 

J n {R) = -2n J dATP(AT, R)[AT\ n ~ 1 H(AT, R) , 

(11) 

where P(AT, R) is the probability density function (pdf ) 
to find a temperature difference AT across a separation 
R. It was proposed by Kraichnan || that the conditional 
average H(AT,R) exhibits in his model a very simple 
functional dependence on AT and on f?, i.e. 



if (AT, R) = 



-J 2 AT{R) 
iS 2 (R) 



(12) 



If this were true, J n {R) could immediately be written in 
terms of structure functions, 



MR) = 



nfoSnjR) 
2S 2 {R) 



(13) 



This again closes the theory upon itself, allowing one 
to proceed. In the Kraichnan model, D n can also be 
written in terms of S n (R) [see Eq.(|5|) below] and one 
can therefore find the scaling exponents that character- 
ize the structure functions. Unfortunately there is still 
no derivation of the ansatz (|l2|). There are indications 
that it is obeyed; numerical simulations of the Kraichnan 
model support it rather convincingly Q. In addition it 
was shown in || , on the basis of experimental data analy- 
sis, that this form of the conditional average is obeyed in 
a context that is much wider than the Kraichnan model. 
Experiments on turbulent advection were analyzed, and 
good agreement with ( |l2| ) was demonstrated. The aim of 
this paper is to present further theoretical and data anal- 
ysis in this direction. We want to understand what can 



be said on conditional averages in terms of fundamen- 
tal theory, and how to intelligently analyze experimental 
data to probe these important quantities. 

It should be pointed out that although we focus in 
this paper on turbulent advection, similar considerations 
are important also in the context of Navier-Stokes turbu- 
lence. Also in that case the two fundamental strategies 
to develop a statistical theory are open to us. The sec- 
ond strategy is even more tempting in that context. In 
the first approach one gets objects depending on one co- 
ordinate, but a hierarchy of equations relating different 
orders (in powers of the velocity field) correlation func- 
tions. The second strategy gives a theory involving many 
coordinates, but in which we can also neglect the viscous 
term, obtaining homogeneous equations T> n = that in- 
volve only one order of correlation functions. We are not 
oing to explore this issue further and refer the reader to 
1| for more details. 

The paper is organized as follows: In Secion II we 
present theoretical considerations that relate conditional 
averages with the probability density function. The fu- 
sion rules are employed to develop a general representa- 
tion of the conditional average ( ^0| ) in terms of the pdf 
of AT(R). It is shown that in general if (AT, R) can 
be written as an expansion in non-integer powers of AT, 
with the first term being linear, and with dimensionless 
coefficients that are denoted <zo,ai .. ., see Eqs.(p2|),(f23"|). 
In Section III we analyze experimental data of passive 
and active scalar advection, with the aim of understand- 
ing whether the linear term in our expansion is leading, 
and whether the rest of the series is unnecessary. We of- 
fer conclusion in Section IV: it turns out that for passive 
scalar advection the linear fits are excellent, whereas in 
the case of active convection the linear form appears to 
fit the data extremely well for high values of the Rayleigh 
numbers Ra. For lower values of Ra there are significant 
nonlinear contributions, and we show that the proposed 
method of data analysis offers excellent fits to the data. 



II. CONDITIONAL AVERAGES AND THE 
RELATIONS BETWEEN THEM 



A. The conditional average of the dissipation field 

In addition to the conditional average (jl^) we will con- 
isder the average of the scalar dissipation field k|VT| 2 
conditioned on AT: 



G(AT, R) 



VT| : 



AT{r,r'] 



(14) 



This conditional average appears naturally in the analy- 
sis of J n (R). For space homogeneous statistics one can 
move one gradient around in the definition (|5j) and get 
for R in the inertial range, 



MR) ~ -2n(n - l)K(\VT(r)\ 2 [AT(r,r')}' 



(15) 
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Accordingly, we can write a second equation in terms of 
the probability density function 

J n (R) = -2n(n- 1) / dATP(AT, R) [AT] n ~ 2 G(AT, R) 



(16) 

By equating @ with (Jll|) wc find an infinite set of inte- 
gral constraints on the conditional averages. This implies 
that the two conditional averages, H and G, must be uni- 
versally related in order to satisfy these constraints for 
any value of n. The required relationship involves the 
pdf P(AT, R) and has the following form: 

H(AT, R)P(AT, R) = —JL- [ G ( Ar , R)P(AT, R)] , 



(17) 



as can be checked by direct substitution. A formula of 
this form has been discussed before in Q . 



B. Fusion rules and their consequences 

An additional constraint on the conditional averages 
can be obtained using the "fusion rules" that have been 
derived recently. These rules serve to find relationships 
between the two fundamental approaches described in 
the introduction. Specifically, the fusion rules address 
the asymptotic properties oiT n when a group of p points, 
p < n — 1 tend towards a common point Tq (|r$ — r$\ ~ p 
for all i < p), while all the other coordinates remain at a 
larger distance R from Tq (|r,; — t*o| ~ R for i > p, and 
R 3> p). For our particular purposes we need to write 
J n (R) as the result of the following fusion process: 

J n (R) = -Inn lim lim V^„(ri,ri . . .r n ,r') 

n— >r r'.^ra+R 

(18) 

The fusion rules that should be used in such cases were 
displayed in great detail in [EJ in the context of Navier- 
Stokes turbulence. They apply identically also to this 
case. Basically it was shown that all the fusions with- 
out gradients in this case have regular limits, relating T n 
with S n . The fusions with gradients require special care 
of the limit ria = r% — r 2 — > 0. The intermediate result 
is, for R in the inertial range, 



J n (R) ~ —Inn, lim 
ri2— >0 



V 2 ri S 2 (r 12 )S n {R)/S 2 (R) 



(19) 



In evaluating this quantity we interpret the limit r\ 2 — > 
as a limit r\ 2 — > r\. This seems natural for large Peclet 
numbers when 77—5-0. It is important however to stress 
that there is a hidden assumption here. We expect the 
function T 2n to change its analytic behavior as a function 
of 7*12 . This change occurs at the viscous crossover scale 
T). The issue is whether this crossover scale is n and R 



independent. That this is so has beenprcwen for Kraich- 
nan's model of turbulent advection 0, and that this is 
not so has been proven for Navier-Stokes turbulence [EJ. 
We believe that this is more generally true in scalar ad- 
vection due to the linearity of the equation of motion (|^) , 
independently of the statistical properties of the driving 
velocity field. The experimental results analyzed in jq] 
strongly indicate that this is the case in a wide context 
of turbulent scalar fields. With this in mind we write 

J n (R) ~ -2Kn[V 2 ri2 S 2 (r 12 )\ ri2=11 }S n {R)/S 2 {R) . (20) 

Using the fact that the mean of the scalar dissipation 
field, denoted e, is evaluated as e ~ K [^ 2 ,S 2 (p)\ p = v ] , and 
also the fact that in the inertial range J 2 (R) = — 4e, we 
write 



MR) 



nC n J 2 S n (R) 
2S 2 {R) 



(21) 



where C n is an as yet unknown dimensionless coefficient, 
but C 2 = 1. Eq.(fjl|) was suggested for Kraichnan's model 
in jlij and derived in || . We proposed in || that it holds 
in a much wider context and showed experimental data 
in support. 

Having the result (^Tj) we see that the scaling exponent 
of J n {R) is fixed as Q n — (2- One way to understand it 
is to assume that indeed Eq.([l2|) is always valid, and to 
use it in Eq.([TT|) to derive this scaling law. In this case 
Eq.(pH) is recovered with the constraint that the coeffi- 
cients C n are all unity and in particular n-independent. 
We need however to allow for the possibility that ( |T^ ) 
is incorrect, and find alternative representations of the 
conditional average that agree with the scaling law (|2l]). 
This is the subject of the next subsection. 



C. Series expansion of the conditional average 

Let us reject for the time being the possibility that 
H(AT,R) is proportional to AT(R) with a coefficient 
depending only on R. Alternatively, let us consider the 
following model expression for H(AT, R): 



if (AT, R) = 



-J a £{ATP(AT,R)} 



AS 2 {R) P(AT,R) 



(22) 



Here we introduce the dimensionless operator C acting 
on the variable AT(R) as a sum of differential operators: 



00 n 

Eap r 

p=0 



(AT) P 



(23) 



In this representation there is the freedom of a countable 
set of dimensionless coefficients a p . 

From the dimensional point of view H(AT, R) in ( g^ ) 
is of the order of AT, but it has a more complicated func- 
tional dependence on AT and i?, expressed in terms of 



3 



the pdf P(AT, R) with the help of operator C. Comput- 
ing J„(r) with H given by (p2|), one gets Eq.(^ll) as we 
should, but with coefficients C n given by 



n-l 



^=E( n : 1 )(-D p -p 



(24) 



p=0 



Here 



(V) 



(n-l)! 



p\(n — p — 1)! 

are binomial coefficients. We have one obvious con- 
straint, i.e. C2 = ao — a\ = 1. One sees that by an 
appropriate choice of a p , an arbitrary dependence of C n 
on n is possible. 

To exemplify the consequences of this extra freedom 
we will analyze next the implications it has on the scal- 
ing exponents of the Kraichnan model of passive scalar 
convected by an infinitely fast Gaussian velocity field. In 
this case the term D n (R) is known exactly, 



D n (R) = 



V d prf-l + Ch d C fD) 



(25) 



where Cfi is the parameter of the model, < Ch < 2, and 
I? is a dimensional constant. Using the balance equation 
and writing S n (R) cx R^ n one computes 



Cn = W(d~ 2) 2 + 2nd£2C n - \{d - 2). 



(26) 



If all coefficients a,( m >i) are zero, then C n = 1 and we 
have Kraichnan's conjecture for Q n : 



Cn = \y/{d- 2) 2 + 2nd( 2 - ±(d - 2) 



(27) 



which dictates the "square-root" asymptotic behavior 
( n — ► \J ndC,2 /2 in the limit n — + oo. Assume now 
that only ao and a\ are nonzero. From (]24|), C n — 
ao — (n — l)ai or 

C n = 1 - (n - 2)aj , (28) 

and if we use this result in the Kraichnan model we get 



Cn = W(d ~ 2) 2 + 2nd( 2 [l ~{n~ 2)ai] - |(d- 2) 



(29) 



Now the asymptotics of Cn in n are linear: Cn oc n^J— a\ 
for n — > oo. Notice that in this case, a\ has to be neg- 
ative. It is interesting to note that the assumption that 
only the first three coefficients, a , a%, and a 2 are nonzero 
would lead to the conclusion that Cn oc n 3 / 2 , which is not 
allowed in view of the Hoelder inequalities for the sca ling 
exponents. Similarly, one cannot truncate the series ( |23| ) 
at any higher term. Hence, only three possibilities are al- 
lowed for this representation of the conditional average: 



(i) only ao is nonzero and we have Kraichnan's expo- 
nents (^7j); In this case we expect to find a linear 
law 



if (AT, R) 



-J 2 



*S 2 (R) 



AT(R) 



(30) 



(ii) only ao and a\ are nonzero and we have the expo- 
nents (|29]); Note that Eq.(p8|) determines the coef- 
ficients C n in this case, and the nagnitude of the 
dimensionless parameter a\ measures the deviation 
of C n from unity. We will see below that in all high 
Re data a\ seems to be smaller than 10~ 2 . 

(iii) there is an infinite set of non-zero coefficients a p . 

It is interesting to ask whether one can come up with 
an example of infintcly many coefficients a p without vi- 
olating any general requirement about the scaling expo- 
nents. In fact this can be easily done. For example, 
choose a p of the following form 

a P = y^&*[l - exp(-a s )] p + 6 pQ u. - 8 pl v , (31) 

s 

and substitute in Eq.(|24|). The result is 

Cn = ^ + ( n ~ IV 
n-l _ 

+ E( ri p )£&( e-Q '- 1 ) p (32) 

p=0 ' s 

Finally it gives: 

C n = fi+ (ri - \)v + £/3 s exp[-(n - l)a s ] , (33) 



wich satisfies the constraint 

C 2 = /x + v + 2J/3 S exp(-a s ) = 1 

s 

In the limit n 3> max[l/a s ] we find 



(34) 



Cm = W( d - 2 ) 2 + 2ndC 2 [M + (n - - §(d - 2) . 

(35) 

This form has again an asymptotic linear dependence of 
C„ on n, but for intermediate values of n these exponents 
differ significantly from (p9|). We do not ascribe partic- 
ular importance to this result, and exhibit only to show 
that to satisfy the consequences of the fusion rules, we 
have in general considerable freedom in the functional 
dependence of Cn on n. 

It is important to understand now that the series $22\) 
is actually an expansion in non-integer powers of AT. As 
such, it is fundamentally different from the series pro- 
posed in ref. Q] which is in integer powers. The non- 
integer powers are dictated by the functional form of the 
pdf P(AT, i?), which in general is non-analytic. In or- 
der to see this clearly, we consider for example a form 
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P(AT, R) that has been found |J to fit very well the ex- 
perimental data for turbulent temperature fluctuations. 
For different separations R, the pdf is described by the 
following stretched-exponential form: 



P{AT,R) = C*(i?)exp \-a(R)\AT\^ 



(36) 



Substitution of this into Eq. ( p2[ ) gives a series in nan 
integer powers of AT which originate from the differenti- 
ation of (AT) 13 . Any attempt to re-expand the series in 
(AT) m P for non-integer f3 in integer powers of AT leads 
unavoidably to a series with zero radius of convergence. 



III. IS THE CONDITIONAL AVERAGE H(AT, R) 
LINEAR OR NONLINEAR IN AT? 

As we already saw, the present state of the theory does 
not allow an ab-initio determination of the functional 
dependence of the conditional average H on AT. Ac- 
cordingly, we turn now to analyzing experimental data 
to shed light on this issue. As explained, the condi- 
tional average if (AT, R) is linear in AT if and only if 
all the coefficients except ao are zero (i.e. the possi- 
bility (i) discussed in last section). For the other two 
cases, H(AT, R) is a nonlinear function of AT. In ear- 
lier work ||, we found that if (AT, R) is close to a linear 
function of AT which implies that a p ,p 7^ are small 
compared to ao- To make more quantitative statements 
we will perform further analysis of the experimental data 
under the assumption that ao and a\ are nonzero. We 
will see that this form fits the data extremely well, and 
that the coefficient a\ is always small, and it appears to 
become smaller when the Reynolds number is increasing 
and when the chosen separation goes into the bulk of the 
inertial interval. Taking ao and a\ as the only nonzero 
coefficients we find 



H(AT, R) 
-J 2 AT 
4S 2 (R) 



(a + 2ai)+aiAT 



d[logP(AT, R)} \ 
dAT J 



(37) 



Thus, the coefficient a\ indeed measures how nonlinear 
H is. Using the form HI for the pdf T(AT, R), H can 
be rewritten as 



H(AT, R) 
-J 2 AT 
4S 2 (R) 



(I + 3 ai ) - a ia (R)(3(R)\AT\ 



0{R) 



(38) 



in which ao = 1 + a\ is used. When P(AT, R) is asym- 
metric, a more general form with different a and (3 for 
AT > and AT < has to be used. Note that if we 
measure AT in units of its standard deviation (as we 
do below in the data analysis), then the combination of 
parameters aia(R)(3(R) gives a direct measure of the im- 
portance of the nonlinear correction in this equation for 
AT = 1. 



A. Linear Fits 



We begin the discussion of experimental data by 
demonstrating that in the case of passive scalar advec- 
tion the linear form of the conditional average H(AT, R) 
is observed to high precision. Firstly we examine the the- 
oretical prediction (f2l|). The results show that to a good 
accuracy C n ~ I for all n and R ||. 

FIG. 1. A 

plotoflog|J 2 „(.R)/(2nK)] vs. log^n)' 1 J 2 S 2n (R)/S 2 (R)\ for 
n — 2 (squares), 3 (triangles), 4 (diamonds), 5 (stars), and 
6 (circles) and R in the inertial range. The data are taken 
from Yale [10]. The line (slope 1 and intercept 0) is not a fit 
but is the theoretical expectation (21) with C„ = 1. 

We use temperature data measured in the wake of a 
heated cylinder jio). Air of speed 5 m/s flowed past a 
heated cylinder of diameter 19 mm (Reynolds number 
= 9.5 x 10 4 ). The temperature was measured at a fixed 
point downstream of the cylinder on the wake centcrlinc. 
The cylinder was heated so slightly that the buoyancy 
term was unimportant and the temperature acted as a 
passive scalar. Temperature was measured as a function 
of time, and we use here the standard Taylor hypothe- 
sis that surrogates time derivatives for space derivatives. 
In doing so we made sure that the viscous scales are 
properly resolved in this data set. In Fig. 1 we display 
J 2n {R)/(2nK) as a function of (2k)" 1 J 2 S 2n (R) / S 2 (R) for 
n varying from 2 to 6, and for various R values in the iner- 
tial range. We see that all the points fall on a line whose 
slope is unity to high accuracy, and whose intercept (in 
log-log plot) is very closely zero. As was pointed out in 
H this good agreement is a confirmation of the validity 
of the fusion rules. It should be stressed that individual 
tests at various values of n as a function of R corrobo- 



rate the same conclusion, i.e. Eq.(21) is supported by the 



experimental data with C 2n being near unity. The most 
sensitive test of the alleged constancy of the coefficients 
C 2n is obtained by dividing J 2n {R) by nJ 2 S 2n {R) / ' S 2 (R) 
for all the available values n and R. The result of such a 
test is shown in Fig. 2. 

FIG. 2. A detailed test of the coefficient C'2 n , see text for 
details. The symbols are the same as in Fig.l. The small 
systematic decrease of C2n with n may be due to insufficient 
accuracy at the tails of the probability density function which 
becomes more important at large values of n. 



FIG. 3. The conditional average H(AT,R) measured from 
the Yale data [10] normalized by the measured value of 
— J2 / ^S2{R) as a function of AT(R) for three different values 
of R measured in units of the sampling time. The differ- 
ent -R values are designated by triangles (R — 16), squares 
(7? = 128) and circles (R = 1024) respectively. 
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We see that all the measured values of Cm are concen- 
trated within the interval (0.75, 1) for separations within 
the inertial interval. Considering the fact that the quan- 
tities themselves vary in this region over five orders of 
magnitude, we interpret this as a good indication for the 
independence of Cm of R and n. The R independence is 
very clear, and is a direct test of the fusion rules. The 
weak n dependence seems to indicate that Cm decreases 
slightly with n; this may arise from the limited accu- 
racy of the data. We are reluctant to make a strong 
claim about the accuracy of 10'th or 12'th order struc- 
ture functions. If wc accept for now the evidence that the 
coefficients C2n in ( pl| ) are n-independent, it must also 
imply that the conditional average H(AT, R) is linear in 
AT. 

In Fig. 3 we present results from the same data set 
that was used above. We show the conditional average 
H(AT,R) as a function of AT(R) for various values of 
R. The line passing through the data points is not a fit, 
but rather the line required by Eq.(p0[). We note that 
points belonging to different values of R fall on the same 
line, indicating that indeed the conditional average H is 
a function of AT(R) times a function of R, and that we 
identified correctly the function of R as — Jz/^SziR). 



B. Nonlinear Fits 

As explained, the linearity of the conditional average 
of H in ST, [Eq.([30])] was not derived from first princi- 
ples. We therefore proceed now to see whether the more 
general form (38) is supported by the data, and whether 
we can bound from above the values of the parameter a\ . 
To estimate a\ from experimental data, we first estimate 



a(R) and (3(R) from the pdf's evaluated from data using 
© then plot 2H(AT, R)S 2 {R) /[J 2 AT] versus |AT|^ fl ). 
The intercept is given by 1 + 3ai. To study how well fl3S| ) 
can represent the data, we subsitute the estimated value 
of <zi into (^) and compare it with the experimental 
data. Since the passive scalar data shown in Fig. 3 agree 
so well with the linear ansatz, we discuss first a case that 
offers a more stringent test of the form ([}8]). To this 
aim we consider data taken from convective turbulence. 
In this case the temperature is an active rather than a 
passive scalar. The data are taken from the well docu- 
mented Chicago experiment Jli"| , |l2| . The experiment was 
performed in a cylindrical box of helium gas heated from 
below, and the Rayleigh number (Ra) can be as high as 
10 15 . The box has a diameter of 20 cm and a height of 
40 cm. The temperature at the center of the box was 
measured as a function of time, and we use the same 
Taylor-hypothesis to surrogate time for the spatial coor- 
dinate. Fig. 4 display the conditional average H(AT, R) 
computed from these data for three different values of 
the Rayleigh number Ra with R measured in units of the 
sampling time. We see that ([38]) is always a good form 
for describing the experimental data. It is interesting to 
examine how the nonlinearity in the conditional average 
depends on Ra and on the value of R. In Table 1 we 
present a compilation of the best fits for the parameters 
for a range of values of Ra and for a range of values of R. 
The results appear to support the following conclusions: 

1. the value of <xi generally decreases when Ra in- 
creases; 

2. the value of a\ is smaller and remains approxi- 
mately constant when the separation R is deep in- 
side the inertial range. 



FIG. 4. The conditional average H(AT, R) as a function of AT for the turbulent convection data. Shown are representative 
fits of the formula (38) at four values of Rayleigh number (Ra) with the separation R measured in units of sampling time using 
Taylor hypothesis. One sees that the linear fit becomes better as Ra increases, and see Table 1 for a quantitative confirmation 



A good way to have a quick estimate of the impor- 
tance of the nonliear term [cf Eq.(|38|)] is to measure a\a(3 
which is the coefficient of the nonlinear term. We see 
that this coefficient decreases significantly when we go 
from Ra= 6.0 x 10 8 to Ra= 5.8 x 10 14 , becoming about 
0.01 in the middle of the inertial range. 

Next we show similar detailed calculations for the pas- 
sive scalar data of Fig. 3. There is a slight complication 
since the pdf's P(AT, R) are not symmeteric in AT. Ac- 
cordingly we need to fit separately the left and right 
branches of the distribution functions, and find the pa- 



rameters aq_, a_, /3+ and /3_, together with the appro- 
priate values of (qi)+ and (qi)-- After doing all this we 
show the fits in Figs. 8-10. It appears to the eye that the 
quality of the fits is not significantly improved compared 
to the linear fit. To see this more clearly we present in 
Table 2 the values of all the parameters involved in the 
fit. It is seen that the values of the parameters a\ ± are 
close to zero, or more precisely a± ± = ± 0.02. The co- 
efficient that measures the importance of the nonlinear 
correction, i.e a\aj3 is of the order of 10~ 4 for all the 
separations in the bulk of the inertial range. 



G 



TABLE I. Fitted parameters for the turbulent convection data, with the Rayleigh number spanning the range 6 x 10 s to 
5.8 x 10 14 . The parameters a and /3 characterize the probability density function P(AT, R), and a\ is the coefficient of the first 
nonlinear contribution to the expansion of the conditional average H(AT, R) in AT. The separation R is measured in units of 
the sampling time. The value of the parameter ai measures the deviation of the coefficients C n from unity, cf. Eq.(28). The 
combination aia/3 measures the deviation of the conditional average from the linear fit at AT = 1, cf. Eq.(38). 



Ra 


R 


8 


16 


32 


64 


128 


256 


512 


1024 


2048 




P 


0.54 


0.64 


0.71 


0.95 


1.06 


1.31 


1.74 


1.84 


— 


6 x 10 


lOOai 


-40.3 


-22.8 


-16.5 


-12.9 


-8 


-17.6 


-7.3 


-16.7 






100aiQ/3 


- 116 


- 57 


- 37.2 


- 24.9 


- 14.3 


- 14.3 


- 4.4 


- 10.0 






P 


0.52 


0.67 


0.85 


1.05 


1.19 


1.41 


1.81 


2.00 


— 


4 x 10 


lOOai 


-24.8 


-14.6 


-12.8 


-5.4 


-7.2 


-6.3 


-3.1 


-16.4 






100aia/3 


-68 


-33.4 


—22.4 


-8.7 


-9.0 


-8.1 


-2.6 


-10.0 






P 




0.61 


0.68 


0.80 


0.91 


1.15 


1.38 


1.55 


1.68 


7.3 x 10 10 


lOOai 




-9.1 


-7.3 


-7.1 


-1.1 


-16.3 


-12.1 


-14.8 


-5.7 




100aiQ/3 




-24.5 


-17.6 


-12.1 


-3.3 


-15.5 


-9.5 


-11.1 


-3.0 




P 




0.60 


0.65 


0.72 


0.84 


0.96 


1.19 


1.50 


1.43 


6 x 10 11 


lOOai 




-6.9 


-2.8 


-14.0 


-12.9 


-10.9 


-13.7 


-4.4 


-14.2 




100aia/3 




-18.8 


-9.1 


-27.1 


-20.6 


-71.0 


-15.7 


-4.7 


- 11.6 




P 




0.58 


0.64 


0.72 


0.83 


0.94 


1.00 


1.25 


1.45 


6.7 x 10 12 


lOOai 




-7.1 


-20.4 


-17.5 


-8.2 


-6.7 


-3.9 


-7.1 


-8.1 




100aiQ/3 




-22.9 


-48.1 


-34.9 


-14.4 


-12.4 


-9.6 


-9.5 


-7.6 




P 




0.61 


0.68 


0.77 


0.86 


0.88 


1.15 


1.26 


1.42 


4.1 x 10 13 


lOOai 




-2.9 


-1.5 


-0.5 


-2.3 


-1.2 


-0.7 


-30.7 


-10.5 




100aia/3 




6.4 


-4.7 


-1.9 


-6.0 


-2.4 


-2.4 


-22.5 


-9.4 




P 




0.59 


0.63 


0.69 


0.85 


0.93 


1.05 


1.31 


1.56 


5.8 x 10 14 


lOOai 




-4.7 


-2.6 


-4.5 


-5.3 


-0.4 


-2.5 


-8.7 


-3.1 




lOOaiQ/3 




-10.9 


-4.7 


-6.0 


-4.6 


-1.0 


-4.3 


-8.9 


-2.8 



TABLE II. The fit parameters for the nonlinear fits in the passive advection data with Reynolds number (Re) = 9.5 x 10 4 . 
The probability density function P(AT, R) is asymmetric in this case, and we fit separately the left (minus subscript) and right 
(plus subscript) wings. All the parameters that measure the deviation from the linear fits are very small. 



R 


8 


16 


32 


64 


128 


256 


512 


1024 


P- 
P+ 


0.94 
1.15 


1.12 
1.32 


1.28 
1.56 


1.54 
1.63 


1.77 
1.75 


1.85 
1.82 


1.76 
1.87 


1.86 
1.82 


10 J (ai)_ 
10 3 (ai)+ 


21 
18 


67 
14 


11 

-1 


12 
19 


-3 
-6 


9 
6 


10 
0.6 


-5 
12 


10 3 (aia/3)_ 
10 3 (aia/3) + 


10 
1.8 


16 
-0.7 


3 

-1.7 


1 

0.2 


-0.4 
-0.3 


0. 
-0.1 


0.2 

-0.2 


-0.4 
0.4 
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FIG. 5. The conditional average H(AT, R) as a function of AT for the passive scalar data at three values of the separation 
R measured in units of sampling time. The nonlinear fits are indistinguishable from the linear ones in the bulk of the inertial 
range, and see Table 2 for a quantitative confirmation 



IV. CONCLUSIONS 

We have presented a theoretical analysis of the re- 
lation between the two conditional averages H(AT, R) 
and G(AT, R) and the probability distribution function 
P(AT,R). The general relation is given by Eq.(17). 
From this relation it follows that if one of these aver- 
ages factorizes to a function of AT times a function of R, 
the other cannot factorize as long as the the distribution 
function does not factorize. The latter cannot factorize 
in any multiscaling statistics. Next we have presented ev- 
idence that the conditional average H(AT, R) does fac- 
torize into a function of R times AT. This appears to 
be the case for both passive and active scalars when the 
Reynolds number is sufficiently large and when R is in 
the bulk of the inertial range. The fusion rules which 
are believed to hold in a wide context furnish a predic- 
tion about the function of R that precedes AT in the 
conditional average, cf. Eq.(f30|). The data support the 
prediction of the fusion rules to very high accuracy. We 
do not have at present a theoretical ab-initio derivation 
of the linear dependence that seems to be supported by 
the data. In view of the importance of this law for the 
study of the balance equation D n (R) — J n (R) it seems 
to us that such a derivation is highly desirable. 
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